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Abstract 

The well-studied local postage stamp problem (LPSP) is the following: given a 
positive integer k, a set of positive integers 1 = ai < < ■ ■ ■ < ak and an integer 

h > 1, what is the smallest positive integer which cannot be represented as a linear 
combination Yli<i<k where i<i<k ^ an d eac h x i is a non-negative integer? 
In this note we prove that LPSP is NP-hard under Turing reductions, but can be solved 
in polynomial time if k is fixed. 


1 Introduction 

The local postage-stamp problem, or LPSP for short, can be informally defined as follows. 
One is given a supply of stamps of k different denominations, 1 = a\ < a 2 < • • • < a*,, and 
an envelope that has room for at most h different stamps. What is the smallest amount 
of postage Nh(a±, ct 2 ,... , a*,) that cannot fit on the envelope? For example, if the available 
denominations are let, 4$, 7(f, and 8<t, and the envelope has room for 3 stamps, then all 
amounts of postage < 24cj: can be provided but 25c cannot. Hence i\r 3 (l, 4, 7, 8) = 25. 

A more formal statement of LPSP is given in the abstract. 

In this note we consider the computational complexity of LPSP. If N = Yh<i<k x i a i, we 
call (xi,X 2 , ■ ■ ■ ,Xk) a representation for N and J2i<i<k x i ^ ie we igiif of the representation. 
If further J2i<i<k x ‘ minimum among all representations for N, we call (x±,X 2 , ■ ■ ■ ,Xk) a 
minimum-weight representation for N. If the denominations a* and bound h are given in 
unary, then a simple dynamic programming algorithm can determine the minimum-weight 
representation for all integers N < hak + 1 in polynomial-time, and hence we can compute 
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Nh(ai, a 2 , ■ ■ ■ , a k ) in polynomial time. We therefore assume for the rest of this paper that 
all inputs are provided in binary. 

LPSP was apparently introduced by Rohrbach |[7]. [8| in 1937, and since then dozens of 
papers have been written about it and a variant, the global postage-stamp problem (GPSP); 
see Guy [[|, pp. 123-127]) for a brief survey. Despite this, no general results on the com¬ 
putational complexity of the problem seem to be known up to now; for example, Alter 
and Barnett [|IJ asked if N h (a i,... , a*,) “can be expressed by a simple formula”. Selmer |J 
discussed efficient algorithms for the case where k < 3. 

In the next section we prove that LPSP is NP-hard under Turing reductions, and in 
Section 3 we give a polynomial-time algorithm for LPSP when k is fixed. 


2 LPSP is NP-hard 


We prove that LPSP is NP-hard by reducing from a related problem, the Frobenius problem 
(see, for example, Guy [||, pp. 113-114]). In the Frobenius problem, we are given an integer 
k > 1 and k positive integers op, a 2 , ... ,a k with gcd(ai, a 2 ,... , a*,) = 1 , and we are asked to 
compute g(ai,a 2 , ... ,a k ), the largest integer which cannot be expressed as a non-negative 
integer linear combination La <i<k x i a i- The Frobenius problem is well-studied, but it was 
only fairly recently that it was proved NP-hard (under Turing reductions) by Ramfrez- 
Alfonsfn ||. 

Before we give the reduction, we need a technical lemma. 


Lemma 1 Let 1 


and 


Gp < a 2 < • • • < a k . Define 

Gp_|_ 1 

K= 2 ^ 




h\ — h o + 


(ho + l)afe-i 





^k —1 


Then 





(a) 

Nho(ai, a 2 ,.. 




(b) 

(^15 d 2 , 

.. , a k ) > (i + 1 )a k for all i > 0; 


(c) 

Nh(a±, a 2 ,... 

, a k ) > (h + l)a fc _i - 

- a k for all h > hi. 


(d) 

Nh+i(ai, a 2 ,. 

• • , a k ) = Nh(ai, a 2 , 

.. , a k ) + a k for all h > hi. 


(e) 

There exists 

a constant c > — 1 such that ha k — A^cp, a 2 ,... , a k ) 

= c for all h > hi. 

(f) 

Ifh> hi, then N h (ai, a 2 ,... , a k ) 

h&k ^(d/c d/c_i, cik ... 

, a*, —ai, a k ), where 


g is the Frobenius number. 




Remark. Parts (a)-(f) can be essentially found more or less verbatim in the paper of Selmer 
the only difference in our presentation is that we explicitly compute the bounds h 0 , hi. 
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Proof, (a) Consider finding a representation n = Yhi<i<k x i a i f° r an integer n, 0 < n < 
a k , using the greedy algorithm. We use at most [a k /a k -i\ copies of a k - 1 , then at most 
[ak-i/cLk- 2 \ copies of a k ~ 2 , etc. The choice of h 0 thus allows us to form the greedy repre¬ 
sentation of all such n. 

(b) We prove this by induction on i. For i = 0 the result is just part (a). Otherwise, 
suppose N ho+i (ai,a 2 , ■ ■ ■ ,a k ) > (i + l)a k . Then every integer n, 0 < n < (i + 1 )ak has a 
representation of weight < h 0 + i and by adding a single copy of a k , we can represent every 
integer m, 0 < m < (i + 2 )ak with weight < h 0 + i + 1. 

(c) Set i = j+ ’ 3 - °> and a PP^ ( b )’ 

(d) The numbers that have representations of weight < h + 1 can be divided into two 
not necessarily disjoint subsets: 

51 = {n : some representation of n of weight < h + 1 has Xf. > 0} 

5 2 = {n : some representation of n of weight < h + 1 has Xk = 0}. 

Now every element of Si can be written as ak + t, where t has a representation of weight 
< h. It follows that 

{ak, ak + 1, ak + 2,... , W(a'i ; • • • , <hc) + — 1} ^ Si, 

but Nh(a±, h 2 , ■ ■ ■ , ak) + ^ S\. On the other hand, the numbers in S '2 have representations 

of weight < h+ l using just the numbers {ai, a 2 , ■ ■ ■ , ak- 1 }, and so the largest clement of S 2 
is < (h+l)ak-i- Furthermore, by (a) and the fact that h\ > h 0 , we have {0,1,... , ak} C S' 2 . 
It follows that provided 

Nh(a\, «2, ■ ■ ■ , Q’k) + ak > (h + l)ak-i, (1) 

we have N h+ i(ai, a 2 ,... , a k ) = W(«i, a 2 ,... , a k ) + a k - But (Jl|) follows from (c). 

(e) Using (d), a simple induction gives N hl+i (ai, a 2 , ■ - - , a k ) = N hl (ai, a 2 ,... , a k ) + ia k 
for all i > 0. Then (hi + i)a k — N hl+i (ai, a 2 ,... , a k ) = hia k — N hl (ai, a 2 ,... , a k ) for all 
i > 0; so we may take c = h\a k — Nhi(ai, a 2 ,... , a k ). Since N hl (ai, a 2 ,... , a k ) < hia k + 1, 
it follows that c > — 1. 

(f) Suppose h > hi. Then by (e) we have ha k — Nh(ai,a 2 , ■ ■ ■ ,a k ) is independent of h. 
Now ha k — t has a representation of weight h' < h iff 

t ha k ^ ^ x^i 

1 <i<k 

= ^2 x i(a k - a^ J + (h - h')a k , 

\l<i<fc—1 / 

i.e., if t has a representation of any weight using the basis a k — a k - 1 , a k — ak - 2 ,... , a k — a 1 , a k , 
since we can choose h to be arbitrarily large. But the largest t with no representation in the 
basis a k — a k ~ 1 , a k — a k ~ 2 , • • • , a k — ai, a k is just the Frobenius number 

d(O i k Ufc—1? Ufc Ufc— 2, • ■ ■ , «1, Ufc)- 

(Since a { = 1, we have gcd(a fc — a k -i,a k — a k - 2 , ■ ■ ■ ,a k — ai,a k ) = 1.) It follows that 
ha k - N h (ai, a 2 , ■ ■ ■ , a k ) = g(a k - a k - 1 , a k - a fc _ 2 ,... , a k - a u a k ). ■ 
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Theorem 2 Given positive integers b\ < b 2 < ■ ■ ■ < b k with gcd(6i, b 2 , • • • , b k ) = 1, 
we can determine in polynomial time integers h,ai = 1, a 2 ,... , a k , a k +i, flfc +2 srtch that 
g{bi, b 2 , • • • j = ha ^+2 — Nh(ai, a 2 ,... , a*,, a^.+i, a*, +2 ). 

Proof. By a theorem of Brauer || Corollary to Thm. 1], we know that g{b\,b 2 ,... ,b k ) < 

Mi- Define 6 fc+ i = Mi and M 2 = Mi +1- Then clearly Mi; b 2 ,... , 6 fc ) = Mi; b 2 , ■ ■ • , h, h+i, h+ 2 )- 
Now we have, by Lemma |l] (f), that 

hiafc+2 — Mi (Mi-2 — bk+ 1, bk+2 ~ bk, ■ ■ ■ , bk+ 2 ~ b\, bk+2 ) = &2, ■ • • , bk, bk+i, bk+ 2 ), 


where 


and 


hi — ho + 


(ho + l)cife+i 
O j k-\-2 ^ k -\-1 


h 0 = 


E 

l<i<fc+2 


®i+l 

Qii 


Note that bk +2 ~ ^fc+i = 1- ■ 

Since the Frobenins problem reduces to LPSP, and the Frobenius problem is NP-hard, 
so is LPSP. 


3 A polynomial-time algorithm for fixed k 

1 observe that results of Kannan [[fj || provide a polynomial-time algorithm for the local 
postage-stamp problem for any fixed dimension k. Let Q be a given copolyhedron (an 
intersection of a finite number of half-spaces, possibly closed, possibly open) in W +l . Let A 
be an m x n matrix, B be a m x p matrix, and C be a column vector of dimension m, all 
with integer entries. Kannan proved that assertions of the form 

“V?/ G Q/T} 3x G 7L n such that Ax + By < C ” 
can be tested in polynomial time when l + n + p is fixed. (Here Qjl} = {y E W : 
there exists 1 £ Z 1 such that [y,t] G Q.) Thus if k is fixed, by taking n — k, l — p — 1, 
and Q = {[y,y] : 0 < y < M — 1}, we can decide in polynomial time whether Vt, 0 < t < 
M- 1, 3xi,x 2 , ... , x k G Z such that x 1 ,x 2 ,... , x k > 0, = b and Y.i<i< k x i - h - 

Now we use a binary search on M to hnd the largest M for which the statement holds. This 
gives us N h (ai, a 2 ,... ,a k ). Note that Q is the intersection of four half-spaces. Kannan’s 
algorithm is quite complicated and this method is likely not to be useful in practice. 
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